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I. WHY THE SUN SHINES 



The weak nuclear interaction plays crucial role in the formation of stars in our Universe: it 
starts the pp chain in stars of the size of our Sun [ll , which is the source of its energy. 
This chain has three branches, ppl, ppll and ppIII [3]. The ppl chain is 

p + p d + + Ue, (1.1) 

p + p + e' d + Ue, (1.2) 

d + p ^ ^He + 7, (1.3) 

^He + ^He ^ ^He + 2p . (1.4) 



In its turn, the chain ppll is 



^He + ^He + 7 , (1.5) 

''Be + ^ ^Li + , (1.6) 
^Li + j9 ^ 2^Ee, (1.7) 



whereas the ppIII chain is 



^Be + p ^fi + 7, (1.8) 
'B 'Be* + e^ + Ue, (1.9) 
^Be -> 2^He. (1.10) 

Besides, the so called hep reaction takes place 

p + ^He ^He + e+ + i/^ . (1.11) 

In these chains, the following reactions occur, triggered by the weak nuclear interaction, 

(1.12) 

p + p + e^ — i- (i+z/e, (1-13) 

(1.14) 
(1.15) 
(1.16) 

The total neutrino flux from the Sun at the surface of Earth is ~ 6.4 x 10^°/cm^ s. The neu- 
trinos produced in the reaction fll.l6p have a continuous spectrum with the maximum energy 
15 MeV. They have recently been registered in the SNO detector 0-0] via the reactions 

fx + d ^ v'^ + n + p , (1-17) 
z/g + e~ + p + (1-18) 

induced by the weak nuclear interaction, too. The measured total flux of active-flavor 
neutrinos is ~ 5 x 10^/cm^s, whereas the flux of electron neutrinos is ~ 1.7 x 10^/cm^s. 
The neutral current to charged current ratio jl] established unambiguously the presence of 
an active neutrino flavor other than z/g in the observed solar neutrino flux, thus confirming 
definitely the phenomenon of the neutrino oscillations and that the neutrinos possess a finite 
mass. When the data from all solar neutrino experiments is combined withjhe KamLAND 
data 0, one obtains ex2=^^At\l degrees and Am?2=7-59lo:2?x 10"^ eV^ 
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II. WEAK REACTIONS IN LABORATORY 



However, the reactions f ll.12p -f ll.18p cannot be studied experimentally with the desired ac- 
curacy in terrestrial conditions at present. In order to grasp them, one should address other 
weak processes in few-nucleon systems that are feasible in laboratories, such as 



n p + e + Vei 
fi' +p n + Uf,, 
fi~ + d — 7- n + n + Ufj^ . 



(2.1) 
(2.2) 
(2.3) 
(2.4) 
(2.5) 



The one-nucleon weak reactions (12. ip and (12. 2p are now experimentally and theoretically 
well explored. The neutron lifetime is < r >woridav. = 879.9 ± 0.9 s ^ and the singlet 
capture rate for the reaction (12.20 . = 725.0 ± 17.4 s~^, has been measured by the MuCap 
Collaboration at PSI T^. As to reactions (12.40 and ^He(/i~, 1/^)^13 (12. 5p . they have been 
studied in great detail as well. As a result, the half-life of the triton is known with an accuracy 



0.3 %, (/Ti/2)i = (1129.6 ± 3) s |ll|, and the capture rate of the reaction ^He(/i 



5 •^^J.J 



Ao= 1496 ± 4 s" Il3l ] is also known with the same accuracy. The situation with the 



reaction ^H(/i , z/^)nn (12. 3p is less favorable so far. Indeed, the last measurements of the 
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As 
which 



doublet capture rate provided Ai/2= 470 ± 29 s"^ [IJ] and Ai/2= 409 ± 40 s" 
we shall discuss in this talk, the experiment planned by MuSun Collaboration [16 
intends to measure A1/2 with an accuracy of ~ 1.5 %, will help to clarify essentially the 
situation in the theory of reactions triggered by weak nuclear interaction in few-nucleon 
systems. 



III. NUCLEAR MATRIX ELEMENT OF THE WEAK INTERACTION 

In order to describe these semi-leptonic reactions one should know how to calculate the 
matrix element of the weak Hamiltonian B.^ between the initial and final nuclear states, 
\i > and \f >, respectively. Such a matrix element enters cross sections and capture rates. 
One can write generally, 

<f\Hw\^>= / dxe''^-^^[j{0) ■ J^{x)f, - UO)J^^,{x)fl (3.1) 

The weak lepton current is the four-vector, given in the V-A theory of the weak interac- 
tions, e.g., for the muon capture 

j^(0) = m(z7)7^(l + 75)M(^). (3.2) 

The Dirac spinor u{k) corresponds to the initial muon with the momentum k and the spinor 
u{i') corresponds to the final muon neutrino with the momentum i7. 

Looking at the matrix element (13.11) we conclude that the problem lies in constructing the 
hadron currents and the nuclear wave functions. In its turn, the isovector one-nucleon weak 
current Jw,fj. is also of the V-A form, 

Jw,M^) = Jv,M^) + JXM)^ (3-3) 
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where the vector part is given by the matrix element of the isovector Lorentz 4-vector current 
operator between the nucleon states, 



9m {qi 
2M 



(3.4) 



and the axial-vector part is analogously, 



gpjqi 

nil 



(3.5) 



M {mi) is the nucleon (lepton) mass and the 4-momentum 
— p^, where p'^ (p^) is the 4-momentum of the final (initial) 



Here a is the isospin index, 
transfer is given by gi^ = p'^ 
nucleon. 

Theoretically, all form factors entering the weak one-nucleon current were well under 
stood and experimentally settled but one It was the induced pseudoscalar gi 

that resisted for quite a long time. Only recently, its value has been fixed in the al- 



ready mentioned MuCap experiment, gp^(g^ = 0.88 m^) = 7.3 ± 1.1 [13, 2^\, which is 



in excellent agreement with the PCAC and chiral perturbation theory (xPT) prediction. 



(g2 = o.88m2) = 8.2 ± 0.2 [20|, 121 



IV. CHIRAL SYMMETRY, SOFT PIONS 

The essence of the problem with the study of the above quoted nuclear reactions is that the 
quantum chromodynamics (QCD) is non-perturbative at low energies. The way of handling 
this obstacle has already been outlined 50 years ago by realizing that the description of 
the electro-weak interaction with a nuclear system, containing nucleons and pions, should 



be based on the spontaneously broken global chiral symmetry SU(2)iX SU(2)j:j |22|, |23 



reflected in the QCD Lagrangian [2^, |25|. This fundamental concept was realized after 
studying the commutation relations of the charges corresponding to the lepton currents for 
the system of electrons and muons with zero masses. In this case, the lepton charges satisfy 
the commutation relations of the group SU(2) x SU(2) (26l |. 
The lepton charges are defined as 

Qi = J (4.1) 
Q51 = J ^^l^l^d^r, (4.2) 

where r is the lepton isospin. 

From Eqs. (14.11) and (14. 2 p one obtains 

[Qii, Qji] = iUjkQki, (4.3) 

[Q5il,Q5jl] = i^ijkQkl, (4.4) 
[Qii , Q5jl] = i^ijkQbkl- (4.5) 

Gell-Mann made an assumption jl^] that the universality of the weak interactions for the 
leptons and hadrons declares itself in that the charges connected with the hadron currents 
satisfy the same simultaneous commutation relations as the lepton currents. 
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Besides, the vector part of the weak nucleon current should satisfy the Conserved Vector 
Current (CVC) hypothesis, that permits the identification of the weak vector current with 
the isovector part of the electromagnetic current. In its turn the weak axial nucleon current 
should satisfy the Partially Conserved Axial Current (PCAC) hypothesis, 

qlJl,{ql) = ^^mlAUql)M:. (4.6) 

Here f-,, is the pion decay constant, rrin is the pion mass, Ap^qf) = l/iml + qf) is the 
pion propagator and is the matrix element of the pion production/absorption amplitude 
between the one-nucleon states. It is seen that the axial current is conserved in the limit of 
zero pion mass. 

This development induced a burst of calculations and powerful low-energy theorems for the 



weak- and electro-production of pions on nucleon at the threshold were established [22|, [23 
Since this concept is correct only for pions with qi = 0, the pions are called soft. Let us 
note that the results based on current algebras are model independent. 
The spontaneously broken global chiral symmetry S\J{2)lX SU(2)/^ predicts the existence 
of three massless particles with the quantum numbers 0~ in the ground state [28]. Since 
this symmetry is broken by the finite mass of the quarks, also these particles acquire finite 
mass, and they can be identified with pions. 

Let us note that Eq. fl3.ip is valid also for the nuclear processes with the electromagnetic 
lepton and nuclear currents. 

V. IMPULSE APPROXIMATION 

Applying Eq. (13.11) in calculations of the processes in nuclei, one first supposed that the 
nuclear current is approximated by the sum of the one-nucleon currents (13.31) . This approx- 
imation is called the Impulse Approximation (lA). Generally, this concept worked well at 
low energies but it was found that, in some cases, it failed to describe the data. First it 
happened in the case of thermal neutron capture by proton, 

n + p ^ d + J , (5.1) 



that the precise experimental value of the cross section, o"*^^^ = 334.2 ± 0.5 mb* [29|, is 
larger by ~ 10 % than the lA cross section, a^^ = 302.5 ± 4.0 mb [30^. 
This reaction is triggered by the space component of the electromagnetic isovector current, 
which is of the order ~ 0{1/M). 

VI. MESON EXCHANGE CURRENTS 

Here for the first time, meson exchange currents (MECs) rescued the situation. The pion 
production amplitude, evaluated in the soft pion limit, provided MECs (see Fig. Ic and 



Fig. Id) that removed about 70 % of the discrepancy |3lH35 . 

Possible presence of the two-nucleon currents follows also from the potential description of 
the two-nucleon system that we consider for the sake of simplicity. In this case, the nuclear 



1 barn (b) = 10" 
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FIG. 1: The possible structure of the two-nucleon current operators; 
Hamiltonian is 

H = T + V, (6.1) 

where T is the kinetic energy and V is the nuclear potential. For the electromagnetic current, 
Jfi{q), the current conservation reads, 

q-J{q) = [H,piq)]. (6.2) 
Writing the current as the sum of the one-nucleon and two-nucleon parts, 

= E U'^^hq^) + U2,q), (6.3) 



one gets 

qi-J{l,i,qi) = [Ti, p{l,i,qi)], 'i = l, 2, (6.4) 
q-J{2,q) = [Ti+Ts, p(2,g)] + ([r,p(l,l,gi)] + (1^2)), (6.5) 

where q = qi + q2- 

It is seen from Eq. (16.51) that this equation cannot be fulfilled if the MECs are absent. 
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VII. CHIRAL LAGRANGIANS, HARD PIONS 



As we have already mentioned, the approach of soft pions is vahd at the threshold energies. 
Next we discuss the concept of chiral Lagrangians allowing to go beyond this restriction. 
The approach of chiral Lagrangians is based on an assumption that any Lagrange theory 
satisfying the requirements 

1. The constructed currents satisfy the fundamental commutation relations known from the 
soft-pion approach. 

2. In the limit of zero pion mass the weak axial current is exactly conserved, 

should reproduce the results of the soft-pion technique jlsj. This is possible thanks to the 
fact that the correct results are obtained already at the level of trees (no loops). The chiral 
Lagrangians of the pion-nucleon system reflect the global chiral s ymm etry. Standardly, they 
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si. This 



are constructed in non-linear realization of the chiral symmetry 
One can also consider the Lagrangians reflecting the local chiral symmetry [3i 
step allowed to extend the nuclear system by the p- and ai mesons, as compensating Yang- 
Mills fields. However, the problem was that the heavy meson masses violated the symmetry. 
Later on, a concept of hidden local symmetry allowed one to avoid this conceptual difficulty 

HSi. 



The nonlinear hard pion chiral Lagrangians [38|, |41|, 142 



served then as starting 
point for constructing the one-boson exchange currents in the tree approximation. We shall 
call this concept as the Tree Approximation Approach (TAA). 

This approach was also apphed to construct the one-boson exchange potentials j43l-47 . 
Besides the potentials of this sort, many phenomenological potentials of various qu ality 



were also constructed 1481-152 
CD-Bonn E 



The precise second generation potentials Nijmegen I |46| . 
and AVI 8 [5^ have ~ 1- 
Let us note that in order to make realistic calculations, both the one-boson MECs and 
potentials should be supplied with the strong form factors by hands. Fortunately, this can 
be done in such a way that the CVC and PCAC hypotheses are still valid. Consequently, 
one can describe the MECs effect consistently: both the MECs and nuclear potentials are 
obtained within the same approach. 

The reactions (fLT71) - (fT:T8|) and were studied within TAA in Refs. These 

reactions are triggered by the space component of the weak current. Its one-nucleon weak 
axial part is of the order ~ 0{1), whereas the space component of the weak axial MECs is 
~ (9(1/M^). This fact makes the calculations of the weak axial MECs effects difficult. 
Here we present for the reaction ^H(/x~, z/^)nn the results, calculated with the first generation 
potentials, for the doublet capture rate A1/2 = 416 ± 7 s~^ [58] and A1/2 = 397.8 - 399.6 
s~^ |5S]. The total effect of the weak MECs was estimated as ~ 30.6 - 33.1 
Recent calculations with the precise potential Nijmegen I provided A1/2 = 416 ± 6 s~^ 



59|. 
6Q|. 



Analogous calculations for the reaction '^He(/i" 

1484 ± 8 s- 



1502 ± 32 s-i 57 and A, 



H provided for the statistical rate Aq 



in a purely phenomenological approach [6I1] . 
Generally, within this approach one can describe well the nuclear phenomena, triggered 



by the electro-weak interaction, up to energies ~ 1 GeV |33M35l. |62|-|68[| . In Fig.|21 we 



present the double differential cross section for the reaction of the backward deuteron 
electro-disintegration 42| . 

e + d^e' + n + p, (7.1) 

where the energetic electrons are scattered backwards, closely to 180 degrees. The detailed 
analysis of the contributions of particular mesons to the cross section has been made in 
Refs. Hlli. 
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FIG. 2: The double differential cross section for the process (jT.ip . 



Prom 1980's, in parallel with the above discussed schemes of the calculations in nuclear 
physics, new approach was worked out, providing more general framework for systematic 
construction of many body currents and potentials. The price for this was explicit absence of 
all degrees of freedom in the Lagrangian but nucleons and pions and applicability to nuclear 
phenomena only in the low energy region. 



VIII. EFFECTIVE FIELD THEORY 



The fundamental step allowing to go beyond the tree approximation was made by Weinberg 
in 1979 ^]. In this work, Weinberg formulated an effective field theory (EFT): if one writes 
down the most general possible Lagrangian, including all terms consistent with assumed 
symmetry principles, and then calculates matrix elements with this Lagrangian to any given 
order of perturbation theory, the result will be the most general possible S-matrix consistent 
with perturbative unitarity, analyticity, cluster decomposition, and the assumed symmetry 
properties. 

So in this way, one can construct dynamical theory, not limited to the tree approximation. 
At low energies, the effective degrees of freedom in nuclear physics are pions and nucleons, 
rather than quarks and gluons. With the heavy mesons and nucleon resonances integrated 



8 



out and the spontaneously broken global chiral symmetry as assumed symmetry, one obtains 
the xPT of the nucleon-pion system. The resulting effective Lagrangian is given by a string 
of terms, dictated by the chiral symmetry, with increasing chiral dimension j72l] . 

This is the QCD Lagrangian of the ttN system at low energies. The Lagrangian of dimension 
one is simple 

611 = ^(^^7m - M + 1^ ^^7^75)* , (8.2) 

where and depend non-linearly on the pion field and external interactions, and 7s are 
the Dirac matrices. Besides, qa is the weak interaction constant. 

At the second order, seven independent terms appear, at the third order, one has 23 inde- 
pendent terms, and at the dimension four, there are 118 independent operators 72|]. Each 
term is multiplied by a (low energy) constant (LEG) that is fixed either in the process of 
the elimination of heavier resonances, or by the data. 



Weinberg also established counting rules [73|, uM allowing one to classify the importance 
of contribution of various diagrams into perturbative expansion of the S-matrix in positive 
powers u of Q/A^, where Q is a quantity characterizing the hadron system (momentum, 
energy or the pion mass) which is small in comparison with the heavy scale ~ 1 GeV. 
As we have already noted above, in order to calculate reliably the capture rates and cross 
sections of reactions one needs to know accurately the nuclear wave functions (potentials) 
and the current operators. 



A. Space component of the weak axial MECs 



The weak axial currents were constructed within the xPT in Refs. [75|, l76[ . At the leading 
(z/=0)- (LO) and the next-to-leading (i^=l) (NLO) orders, the weak nuclear current consists 
of the well-known single nucleon terms. The space component of the weak axial MECs, 
which is of the main interest here because it enters in calculations of observables in all weak 
reactions mentioned above, appears at N^LO {u = 3). We already know that the single 
nucleon current is known precisely. On the other hand, in the space component of the weak 
axial MECs appears one LEC, called d^. This parameter manifests itself in the two-nucleon 
contact vertex with the weak axial current (see Fig.[3]). Besides, it is present also in a contact 
term part of the three-nucleon force constructed within the xPT. It follows that if this LEC 
will be fixed in one of the few-nucleon processes feasible in the laboratory, one can make 
model-independent predictions for other weak processes triggered by this component of the 
weak axial current. 



B. The hybrid calculations 

Up to present only few calculations fulfil the requirement of consistency and so called 'hy- 
brid' approach is applied instead: the current operator is constructed within the xPT as 
outlined above, but the wave functions are generated either from the one-boson-exchange- 
or phenomenological potentials of the TAA approach. Besides, in calculating the observ- 
ables for the three-nucleon processes, also the three-nucleon forces Tucson- Melbourn[77i] and 
Urbana IX [78[ were used. 
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a 



b 



FIG. 3: The general structure of the two~nucleon weak axial operators; a- the long-range operator, 
b- the short-range operator. 



So far almost all calculations aiming to study the weak interaction in few-nucleon systems 
profited from the precise knowledge of the half-life of the triton to extract the LEG d^. In this 
way, in Ref. (79|, this constant was extracted from the Gamow- Teller (GT) matrix element 
for the reaction (12. 4p . and then the spectroscopic factors^ Spp(0)=3.94x (1±0.004) x 10"^^ 
MeV b and S/iep(0) = (8.6±1.3)xl0~'^° keV b were calculated for reactions of the proton- 
proton fusion fll.l2p and the hep reaction (11.141) . respectively. 

Using the same values of d^, in Ref. |80[ the doublet capture rate, Ai/2=386 s^^, for the 
reaction ^H(/i~, z/^)nn was obtained and in Ref. 8l| the cross sections of the reactions. 



(I1.17P and (I1.18p . were calculated. Similarly, in R ef. l82l| . the capture rate Ao= 1499±16 s~^ 
for the reaction ^He(/i~, t'^)^!! was gained. In Ref.[83|, Marcucci and PiaruUi used AV18 NN 
potential to generate the nuclear wave functions for the process ^H(/i~, z/^)nn and AV18 -|- 
Urbana IX NNN potentials for the reaction ^He(/i~, i^„,)^H. The weak current was taken from 



the xPT approach with the potential current of Ref. [84| added. The calculations resulted 
in Ai/2 = 393.1 ± 8 s'^ and Aq = 1488 ± 9 s'^ 

The problem with these calculations is that in the three-nucleon systems, the consistent 
calculations require to know two LEGs, cd and ce |85|. The constant c^ is related to the 
constant d^ as 

^-cd + \mn (c3 + 2c4) + ^ . (8.3) 



d 



R 



K9A 
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Here the LEGs £3 and £4, together with ci, C2 and ce, are obtained either from the vrN- 



or NN scattering [gl 
xPT Lagrangian 



or in the process of elimination of higher resonances from the general 



In the hybrid calculations of the process ^He(yU 



mentioned above, it is not possible 



to fix Cd and ce simultaneously, because ce does not enter them. Besides, the constant cd 
{d^) enters not only the weak axial MEGs, but also the contact and one-pion exchange part 
of the three-nucleon force constructed within the xPT, whereas the constant ce enters the 

^ The spectroscopic factor is defined as S{E) — a{E) E e'^'^''' , where E is the energy, a is the cross section, 
and T] is the Sommerfeld parameter describing the barrier penetrability. 
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three-nucleon contact term ^ 

Let us note also the work [601], where the capture rate A1/2 = 416 ± 6 s^^ for the reaction 
^H(/i~, Vy^mi was calculated within the TAA in various current models and from this rate, 
the model dependence of the value of in hybrid calculations was studied. 



C. Chiral potentials 

The early theory of nuclear force also started from the nucleon-pion system in 1950s but 
failed to describe reasonably well the empirical data. As we have already discussed, only 
later the chiral symmetry was recognized as proper symmetry of the Nature. In 1960s, the 
heavy mesons were discovered and used to construct quite successful models of nuclear force. 
The price for it was the necessity to introduce phenomenological strong form factors into the 
one-boson-exchange potentials and to deal with a not well known scalar-isoscalar a meson. 
With the advent of QCD, the attempts to construct the nuclear forces from the QCD- 
inspired quark models appeared. But only within the xPT it was possible to derive con- 
sistently the nuclear force between the N nucleons from the same Lagrangian and precise 
two-nucleon potentials have been constructed up to N'^LO: the Entem-Machleidt (EM) [s^l 
and Epelbaum-Gloeckle-MeiBner [91!] potentials. In both cases, the entering parameters are 
standardly extracted from the fit to the nucleon-nucleon scattering data and the deuteron 
properties, or some of them are adopted from the analysis of the ttN scattering (911] . 



The three-nucleon force has been derived within the xPT, but only up to N LO |89l . 190 



far. As we already know, this force contains two constants, c^ and c^;, to be determined. 
The first consistent extraction of the constants c^i and ce from the three-nucleon system 



has recently been made by Gazit, Qualioni and Navratil in Ref. j85[, where these constants 
are constrained by simultaneous calculations of the binding energies of the three-nucleon 
systems and of the triton /3 decay. The nuclear wave functions are generated in accurate 



ah initio calculations using both the two-nucleon EM 87|] and three-nucleon N^LO force 
whereas the process (12.41) is calculated with the weak axial MECs derived from the 
same xPT Lagrangian [82;| as the nuclear forces. The resulting values of the two constants 
are restricted in the intervals, -0.3 < c/j < -0.1, -0.220 < c^; < -0.189. 
The reactions of muon capture in deuterium and in ^He have very recently been calculated 
by the Pisa group (MEAL) j93| in the TA- and hybrid approaches and within the xPT 
as well. The resulting values of the capture rates are, A1/2 = (389.7 - 394.3) s~^ and 
Ao = (1471 - 1497) s-^ 

In order to compare our results with those of MEAL, we choose the last row of TABLE VI 



93l | in which we divide the items by 1.024 |2(J]. This number represents the inner radiative 



corrections taken into account by MEAL. In our analogous yPT calculations ^] with the 



same N^LO potential [STJ, we take the value of c/j = -0.2 85j and A = 500 MeV. Besides, 
using the values of the LECs (in unit GeV~^), Ci = -0.81, C2 = 2.80, C3 = -3.20 and C4 = 
5.40 [ST^, we obtain from Eq. ( I8l3|) = 2.400. Let us note that this value of d^ differs 



considerably from d^ = 1.00(9), presented in the last row of TABLE V 9^ 
The results of calculations are given in Table [H As it is seen, the main difference stems from 
the ^So and ^Pi waves. Also the contributions for the wave differ, but they are small. 
Since the calculated MECs effect, AMEC = 15.4 s^^ 941, is similar to the one obtained in 



[60|, l80|, we conclude that the difference comes from the lA calculations in the channels ^Sq 
and ^Pi. 
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TABLE I: Contributions to Ax/2 from partial waves (in s ^). In the last column, all the contribu- 
tions are summed up. RT - our calculations, MEAL - the last column of TABLE VI [93|] divided 
by 1.024. 





iSo 3p^ total 


RT 

MEAL 


258.7 19.6 59.5 69.7 6.2 0.4 414.1 
244.6 19.4 45.3 69.8 4.3 0.9 384.3 



IX. CONCLUSIONS 

We have seen that the constant [ce) is currently extracted from the triton beta decay 
rate. However, it is abundantly clear from our discussion that dealing with the complexity 
of the three-nucleon system can be avoided, if would be determined accurately from the 
muon capture in deuterium. Then one can consistently calculate other two-nucleon weak 
processes, such as proton-proton fusion reaction fll.l2p and both reactions, fll.l7p and fll.lSp . 
of solar neutrinos with the deuterons, so important for the astrophysics. Besides it turns 



out 92| that d^ enters also the capture rate for the reaction vr )'7 + 2?t,, which is 



the best source of information on the neutron-neutron scattering length a^ 
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